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PROGRAM  INFORMATION 


Description  of  program 

CSDOOR  called  X0057  In  the  ^nversationally  lOrlented  teal-Tlme 
programming  System  (CORPS)  library*  is  useful  to  perform  rapid  design  of  metal 
plates  used  to  form  the  sides  and  roofs  of  blast  cells  and  of  metal  used  as 
doors—regular  and  built-up.  The  program  may  be  used  for  any  structural 
materials  for  which  the  material  properties  are  known,  but  steel  is  the  most 
commonly  used.  The  program  may  be  used,  with  limitations,  to  optimize  the 
design  by  finding  the  least-cost  structures  that  satisfy  all  the  design 
constraints . 

Coding  and  data  format 

CSDOOR  is  written  in  FORTRAN  and  is  operational  on  the  following 
systems : 

a.  U.  S.  Army  Engineer  Waterways  Experiment  Station  (WES)  Honeywell 
DPS/1. 

b.  Office  of  Personnel  Management  Honeywell  6000  Series  at  Macon,  Ga. 

£.  Control  Data  Corporation  CYBERNET  System,  CYBER  170/760. 

Data  can  be  input  either  interactively  from  the  user's  termitial  or  from 
a  prepared  data  file  with  line  numbers.  When  data  are  input  from  the 
terminal  during  execution,  the  program  provides  prompting  messages  to  indicate 
the  type  and  amount  of  input  data  to  be  provided.  Output  can  be  obtained  at 
the  terminal,  written  to  a  permanent  file  for  listing  at  the  terminal  at  a 
later  date,  or  directed  to  a  mainframe  line  printer. 

How  To  use  CSDOOR 

A  short  description  of  how  to  access  the  program  on  each  of  the  three 
systems  is  provided  below.  It  is  assumed  that  the  user  knows  how  to  sign  on 
the  appropriate  system  before  trying  to  use  CSDOOR.  In  Che  example  initiation 
of  execution  comnands  below,  all  user  responses  are  underlined,  and  each 
should  be  followed  by  a  carriage  return. 

WES  DPS/4  and  Macon  systems 

After  the  user  has  signed  on  the  system,  Che  two  system  coasands  FORT 
and  NEW  get  the  user  to  the  level  to  execute  the  program.  Next,  the  user 
Issties  the  run  command 

RUM  WESLIB/CORPS/X00S7.R 
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to  Initiate  execution  of  the  program.  The  program  Is  then  run  as  described 
In  this  user's  guide.  The  data  file  should  be  prepared  prior  to  issuing  the 
RUN  command.  An  example  of  initiation  of  execution  is  as  follows,  assuming  a 
data  file  had  previously  been  prepared: 

HIS  SERIES  600  ON  10/17/83  AT  CHANNEL  5647 
USER  ID  -  R0KACASEMP 
PASSWORD  - 
SYSTEM?  FORT  NEW 

0364200 

COEWES  HIS  TIMESHARING  ON  10/17/83  AT  10.328  CHANNEL  2244  TSl 

USER  ID  —ROKAOMP 
PASSWORD  ~ 

USERS  =  043  TSS=145K  MEM-USED=38  SYS-0128K  PRO=2  OOO-WAIT-OOOK 

*F0RT  NEW 
READY 

*RUN  WESLIB/CORPS/X0057.R 


CYBERNET  System 

The  log-on  procedure  is  followed  by  a  call  to  the  CORPS  procedure  file 
OLD , CORPS /UN»CECELB 


to  access  the  CORPS  library.  The  file  name  of  the  program  is  used  in  the 
command 

BEGIN, CORPS, X0057 


to  initiate  execution  of  the  program.  An  example  is: 
83/10/17.  11.07.05.  AA313IA 

EASTERN  CYBERNET  CENTER  SN487  NOS  1. 4/531. 281/18AD 

FAMILY:  KOE 
USER  NAME:  CEROOO 

TERMINAL:  515,  NAMIAF 
RECOVER/CHARGE:  CHARGE, CEROXXX, CEROOO 
$CHARGE , CEROXXX , CEROOO 

07.05.42.  WARNING 

ALL  KOE  USERS  PLZ  TYPE  EXPLAIN, RBF, CONFIG. 
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C>OLD . CWtPS/DN-CECELB 
C>BEGIN. . CORPS. X0057 


How  To  Use  CORPS 

The  CORPS  system  contains  many  other  useful  programs  which  may  be 
catalogued  from  CORPS  by  use  of  the  LIST  command.  The  execute  command  for 
CORPS  on  the  WES  and  Macon  systems  Is: 

RUM  WESLIB/CORPS/CORPS.R 

ENTER  CQMfAMD  (HELP,LIST, BRIEF, MESSAfiE.EXECUTE, OR  STOP) 

*?L1ST 

on  the  Boeing  system,  the  commands  are: 

OLD .CORPS/UN-CECELB 

ENTER  COMMAND  (HELP,L1ST,BR1EF,MESSAGE,EXECUTE,0R  STOP) 

*?LIST 
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PREFACE 


This  user's  guide  documents  a  computer  program  called  CSDOOR  that  can 
be  used  for  optimum  dynamic  design  of  nonlinear  metal  plates  subjected  to 
blast  (pressure-time)  loading.  CSDOOR  is  a  modified  version  of  a  program 
called  SDOOR  that  was  written  by  Mr.  John  M.  Ferrltto,  Civil  Engineering 
Laboratory  (CEL),  U.  S.  Navy  Construction  Battalion  Center,  Port  Hueneme, 
Calif.  The  SDOOR  program  was  modified  to  include  peak  gas  pressures  used  by 
the  Huntsville  Division  (fig.  3-9  described  in  HNDM-1 11 0-1-2)  and  to  allow 
it  to  execute  in  a  time-sharing  mode  with  free-fleld  input.  The  program  is 
useful  to  perform  rapid  design  of  metal  plates  which  function  as  sides  and 
roofs  of  blast  cells  and  chose  which  serve  as  doors — regular  and  built-up. 
Detailed  analysis  and  the  fine  points  of  design  should  be  considered  by 
other  means.  For  example,  design  of  connections  should  be  in  accordance 
with  Huntsville  Division  Report  HNDH-1110-1-2,  "Suppressive  Shields  Structural 
Design  and  Analysis  Handbook." 

Funds  for  modification  of  Che  SDOOR  program  and  preparation  of  this 
user's  guide  were  provided  Co  the  Automatic  Data  Processing  (ADP)  Center, 

U.  S.  Army  Engineer  Waterways  Experiment  Station  (WES),  by  the  Office,  Chief 
of  Engineers  (OCE),  under  the  Computer-Aided  Structural  Engineering  (CASE) 
Project. 

The  program  was  tested  and  reconraended  for  Corps  of  Engineers'  use  by 
Che  CASE  Task  Group  on  Structures  Subject  to  Explosion: 

Mr.  Robert  Wamsley,  Huntsville  Division  (Chairman) 

Mr.  Dennis  Bellet,  Sacramento  District 

Mr.  John  M.  Ferrltto,  U.  S.  Navy  Construction  Battalion  Center 

Mr.  Byron  Foster,  South  Atlantic  Division 

Mr.  William  Gaube,  Omaha  District 

Mr.  William  Hill,  Middle  East  Division 

Dr.  Paul  F.  Mlakar,  WES 

Major  parts  of  this  user's  guide  are  taken  directly  from  Mr.  Ferrltto' s 
original  report  on  SDOOR  (CEL  Technical  Note  TN-1518).  Mr.  Wamsley  and 
Mr.  Paul  K.  Senter,  ADP  Center,  WES,  were  responsible  for  the  program  modifi¬ 
cations  and  wrote  those  parts  of  the  report  pertaining  to  the  modifications. 
Dr.  N.  Radhakrishnan,  Chief,  ADP  Center,  WES,  monitored  the  work,  assisted 
by  Mr.  Senter.  Mr.  Seymour  Schneider  and,  later,  Mr.  George  Matsumura, 
Military  Programs  Directorate,  were  the  OCE  points  of  contact. 
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Commanders  and  Directors  of  WES  during  the  period  of  development  were 
COL  Nelson  P.  Conover,  CE,  and  COL  Tllford  C.  Creel,  CE.  Technical  Director 
was  Mr.  Frederick  R.  Brown. 
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CONVERSION  FACTORS,  U.  S.  CUSTOMARY  (NON-SI)  TO  METRIC  (SI) 
UNITS  OF  MEASUREMENT 


U.  S.  customary  (NON-SI)  units  of  measurement  used  In  this  report  can  be 
converted  to  metric  (SI)  units  as  follows: 


Multiply _  By  _ To  Obtain 


cubic  feet 

0.02831685 

cubic  meters 

Fahrenheit  degrees 

5/9 

Celsius  degrees* 

feet 

0.3048 

meters 

inches 

2.54 

centimeters 

pounds 

(force)  per 

inch 

1.75126850 

newtons  per  centimeter 

pounds 

(force)  per 

square  inch 

6.89475789 

kilopascals 

pounds 

(mass) 

0.45359237 

kilograms 

square 

feet 

0.09290304 

square  meters 

*  To  obtain  Celsius  (C)  temperature  readings  from  Fahrenheit  (F)  readings, 
use  the  following  formula:  C  ■  (5/9)  (F  -  32). 


5 


USER'S  GUIDE:  COMPUTER  PROGRAM  FOR  OPTIMUM  DYNAMIC 
DESIGN  OF  NONLINEAR  METAL  PLATES  UNDER 
BLAST  LOADING  (CSDOOR)* 


PART  I :  INTRODUCTION 
Computer  Program  CSDOOR 

1.  Program  CSDOOR  was  developed  to  perform  rapid  design  of  metal 
plates  used  to  form  the  sides  and  roofs  of  blast  cells  and  of  metal 

plates  used  as  doors.  The  program  includes  provisions  for  the  use  of  plastic 
section  modules  for  built-up  doors,  but  optimization  of  such  doors  should 
not  be  performed  because  the  weight-strength  function  is  not  defined. 
Nonhomogeneous  door  sections  are  not  usually  thought  of  as  plates.  However, 
for  simplification,  the  term  "plate"  will  be  used  throughout  this  guide  to 
denote  these  sections,  whether  they  be  sides,  roofs,  walls,  or  built-up 
doors.  Also,  the  program  may  be  used  for  any  structural  material  for  which 
the  material  properties  are  known;  however,  because  steel  is  more  commonly 
used  in  construction,  the  word  "steel"  will  be  used  in  this  guide. 

Background  of  Original  Computer  Program 
(SPOOR)  Development 

2.  The  Department  of  Defense  (DOD)  has  numerous  facilities  in  which 
various  types  of  explosives  and  munitions  used  by  military  services  are 
produced.  In  most  cases  the  production  of  ammunition  utilizes  assembly  line 
procedures.  Projectiles  pass  through  various  stages  of  preparation:  filling 
with  explosive,  fuzing,  marking,  and  packing.  Hazardous  operations,  such  as 
the  filling  of  the  projectile  case  with  an  explosive  in  a  powder  form  and 
the  compaction  of  the  powder  by  hydraulic  press,  are  accomplished  in  protec¬ 
tive  cells  that  are  Intended  to  confine  the  effects  of  an  accidental  explosion. 


*  CSDOOR  is  designated  X0057  in  the  Conversationally  Oriented  Real-Time 
£rogram-Generating  System  (CORPS)  library.  Three  sheets  entitled 
"Program  Information"  have  been  hand-inserted  inside  the  front  cover 
of  this  report.  They  present  general  information  on  the  program  and 
describe  how  it  can  be  accessed.  If  procedures  used  to  access  this 
and  other  CORPS  library  programs  should  change,  recipients  of  this  re¬ 
port  will  be  furnished  a  revised  version  of  the  "Program  Information. " 
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3.  Most  of  the  existing  production  facilities  were  built  in  the 
1940 's.  With  few  exceptions,  the  manufacturing  technology  and  existing 
equipment  represent  the  state  of  the  art  at  that  time.  The  production 
equipment  was  operated  extensively  during  World  War  II,  again  during  the 
Korean  conflict,  and  recently  during  the  Southeast  Asia  war.  Much  of  this 
equipment  and  many  of  the  structures  housing  it  have  been  operating  beyond 
their  designed  capacities. 

4.  DOD  is  conducting  an  ammunition  plant  modernization  program 
intended  to  greatly  enhance  safety  in  the  production  plants  by  protective 
construction  and  automated  processing  and  to  reduce  the  number  of  personnel 
involved  in  hazardous  operations  (Mendolia  1973). 

5.  In  1969  Technical  Manual  5-1300  (Departments  of  the  Army,  Navy, 
and  Air  Force  1969)  was  published  to  provide  guidance  to  structural  designers 
of  munition  plants.  Specific  objectives  of  the  manual  were  to  establish 
design  procedures  and  construction  techniques  to  (a)  prevent  propagation  of 
explosions  from  one  building  (or  part  of  a  building)  to  another,  (b)  prevent 
mass  detonations,  and  (c)  protect  personnel  and  equipment.  The  manual 
establishes  blast-load  parameters  for  designing  protective  structures  and 
provides  methods  for  calculating  the  dynamic  response  of  steel  plates.  The 
design  method  accounts  for  close-in  effects  of  a  detonation  with  its  associated 
high  pressures  and  nonuniformity  of  loading  on  protective  barriers. 

6.  A  detailed  method  for  assessing  the  degree  of  protection  afforded 
by  a  protective  facility  did  not  exist  prior  to  this  manual's  publication; 
consequently,  the  manual  represents  a  significant  improvement  in  design 
methods.  The  simplifications  made  in  the  development  of  the  design  procedures 
are  presented  in  the  manual.  The  analysis  of  a  structure  using  the  design 
procedure  will  generally  result  in  a  conservative  estimate  of  the  structure's 
capacity;  therefore,  structures  designed  using  these  procedures  will  generally 
be  adequate  for  blast  loads  exceeding  the  assumed  load  conditions. 

7.  Even  with  the  simplifications  presented  in  TM  5-1300,  the  computa¬ 
tional  procedures  are  complex  and  time-consuming.  An  automated  procedure 
was  needed  to  enable  structural  designers  to  perform  rapid  analyses  of  the 
structural  safety  of  blast-resistant  walls  and  doors.  The  design  parameters 
Interact  in  a  complex  way  since  the  procedure  is  both  nonlinear  and  dynamic. 
From  a  design  point  of  view,  an  optimization  procedure  was  needed  to  minimize 
cost  and  maximize  safety  since  blast-resistant  construction  has  been  reported 
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to  cost  three  to  five  times  as  much  as  conventional  construction.  Therefore, 
the  first  objective  was  to  automate  the  analysis  procedures  for  determining 
the  structural  response  of  steel  plates  having  a  bilinear  stiffness  represen¬ 
tation  and  subjected  to  blast  shock  and  gas  pressures.  Plates  are  the  basic 
elements  forming  sidewalls,  roofs,  floors,  and  doors  of  cells  designed  to 
confine  the  effects  of  accidental  explosions.  The  second  objective  was  to 
provide  an  optlmiom  design  procedure  that  will  automatically  produce  a  least- 
cost  design  for  a  given  geometry  as  well  as  material  properties  and  explosive 
weight  for  both  feasible  and  nonfeaslble  starting  points. 

8.  To  meet  these  objectives,  the  U.  S.  Navy  Construction  Battalion 
Center  developed  a  program  SDOOR  that  automated  the  analysis  procedures  and 
optimized  the  design  procedures  (Ferrltto  1977).  However,  additional  modifi¬ 
cations  to  the  program  were  needed  to  Include  peak  gas  pressures,  which  are 
used  by  the  U.  S.  Army  Engineer  Division  Huntsville  In  design,  and  to  allow 
the  program  to  execute  In  a  time-sharing  mode  with  free-fleld  Input.  Personnel 
of  the  U.  S.  Army  Engineer  Waterways  Experiment  Station,  the  Huntsville 
Division,  and  the  U.  S.  Navy  Construction  Battalion  Center  collaborated  to 
modify  SDOOR— hence  CSDOOR. 
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PART  II:  ANALYTICAL  METHODS  AND  LIMITATIONS 


9.  In  general,  the  analytical  methods  used  in  computer  program  CSDOOR 
follow  those  In  TM  5-1300;  consequently,  the  accuracy  of  both  is  the  same. 

The  methods  are  discussed  in  detail  in  TM  5-1300,  TN  1494  (Ferritto  1977)  and 
WES  Instruction  Report  K-81-6,  neither  of  which  will  be  presented  here.  The 
CSDOOR  solution  of  Che  dynamic  response  equation  of  motion  has  been  found  to 
agree  very  closely  with  the  response  chart  of  TM  5-1300.  Additionally,  the 
solution  covers  a  wider  range;  thus  it  is  more  accurate  in  the  areas  not 
defined  by  the  response  chart.  When  the  loading  is  less  than  one  twentieth 
of  the  natural  period,  the  response  is  determined  by  impulse  equilibrium. 

The  basic  dynamic  model  in  TM  5-1300  is  limited  to  the  primary  response  mode 
and  does  not  consider  higher  modes. 

10.  The  blast  impulse  computation  is  restricted  to  a  wall  geometry  in 
which  the  height-to-length  ratio  is  greater  than  0.2.  A  modification  was 
made  by  the  Naval  Surface  Weapons  Center  to  the  original  Picatinny  Arsenal 
Program  to  remove  several  minor  problem  areas,  such  as  the  location  of  the 
charge.  The  blast  impulse  has  all  the  limitations  associated  with  the 
original  Picatinny  programs  that  are  caused  by  limitations  in  the  test  data. 
It  assumes  the  charge  is  an  equivalent  sphere  of  TNT.  Shape  effects, 
explosive  equivalence,  and  explosive  casings  are  considered,  but  only  in  an 
empirical  manner  as  a  result  of  limited  available  data. 
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PART  III:  STRUCTURAL  OPTIMIZATION 


11.  The  cptlmlzatlon  problem  consists  of  finding  the  least-cost 
structure  that  satisfies  all  the  design  constraints;  or,  stated  in  optlmlza- 
tlon  terms,  it  consists  of  finding  X  such  that  M(X)  is  a  minimum  and 

g^(X)  <0  1  -  1,  2,  N 

where 

^  >  vector  of  design  variables 
N  >■  number  of  design  constraints 
g  *  vector  of  design  constraints 
M  *  objective  function 

Specifically  for  this  problem,  the  design  variable  selected  is  thickness  of 
steel  plate  or  section  properties  for  a  nonhomogeneous  section.  The  design 
constraints  are  the  flexural  and  shear  limits.  The  objective  function  is 
the  cost  of  the  steel. 

12.  The  fixed  variables  are 
W  •  explosive  weight 

H  ■  height 

EL  ■  length 

h  »  height  of  explosive  above  floor 

i  *  distance  of  explosive  from  left  side  of  wall 

R  ■  distance  of  explosive  from  wall 
a 

I  *  reflection  code 
f  *  dynamic  yield  stress 
ji  »  ductility 

13.  The  design  parameter  X  is 

X  ■  t  (thickness  of  plate) 

14.  The  constraints  g(X)  are 

6(X)  *  6(9)  ,  maximum  deflection 
t  i  0.05  ,  minimum  tl  <  ?s 
t  £  20  ,  maximum  th i l 

15.  The  methodology  selc  he  unconstrained  minimization 

approach  (Fox  1971  and  Pope  1971  n  problem  is  converted  to  an  uncon¬ 
strained  minimisation  by  constructing  a  function  ^  of  the  general  form 


(K^.r)  »  M(?)  +  P[gj^(i^) . 


For  this  problem  the  Interior  penalty  function  technique  was  selected.  This 
methodology  is  suitable  when  gradients  are  not  available,  and,  because  the 
method  uses  the  feasible  region,  a  usable  solution  always  results.  The 
objective  function  is  augmented  with  a  penalty  term  that  is  small  at  points 
away  from  the  constraints  in  the  feasible  region  but  increases  rapidly  as 
the  constraints  are  approached.  The  form  is  as  follows: 

N 

*(^,r)  »  M(X)  -  ry^  — i— 
j^  8j(X) 

where  M  is  to  be  minimized  over  all  ^  satisfying  g(5)  <  0  ,  j  =  2  , 
...,N  .  Note  that  if  r  is  positive,  then,  since  at  any  interior  point  all 
of  the  terms  in  the  sum  are  negative,  the  effect  is  to  add  a  positive  penalty 
to  M(^)  .  As  the  boundary  is  approached,  some  g(X)  will  approach  zero, 
and  the  penalty  will  increase  rapidly.  The  parameter  r  will  be  made 
successively  smaller  in  order  to  obtain  the  constrained  minimum  of  M  . 

16.  As  stated  above,  the  objective  function  (F)  is  the  cost  of  the 
steel.  Or 

Cost  =  F  =  H*EL*fC 


where  C  is  the  volumetric  cost  of  material. 


N 

•t*  =  F  +  r 

j=l 


1 

gj(^) 


where  r  is  the  penalty  parameter. 

17.  The  program  requires  a  starting  point  in  the  feasible  region 
before  optimization  can  proceed.  The  program  can  automatically  determine 
the  starting  point  by  incrementing  the  design  variables  until  a  feasible 
point  is  reached. 

18.  An  algorithm  which  comprises  the  steps  most  commonly  used  is  as 
follows : 

a.  Given  a  starting  point  X^  ,  satisfying  all  gj (^)  ^  0  » 
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and  an  initial  value  for  r  ,  minimize  if)  to  obtain 

X  .  . 

min 

Qieck  for  convergence  of  to  the  optimum. 

£.  If  the  convergence  criterion  is  not  satisfied,  reduce  r 
by  r  rc  ,  where  c  <  1  . 

Compute  a  new  starting  point  for  the  minimization, 
initialize  the  minimization  algorithm,  and  repeat  from 
step  a. 

The  logic  diagram  for  the  Interior  penalty  functions  technique  is  shown  in 
Figure  1. 

19.  The  minimization  for  ^(^,r)  shown  in  Figure  1  is  accomplished 
by  a  method  developed  by  Powell  using  conjugate  directions  (Fox  1971). 
Powell's  method  can  be  understood  as  follows:  Given  that  the  function  has 
been  minimized  once  in  each  of  the  coordinate  directions  and  then  in  the 
associated  pattern  direction,  discard  one  of  the  coordinate  directions  in 
favor  of  the  pattern  direction  for  inclusion  in  the  next  m  minimizations, 
since  this  is  likely  to  be  a  better  direction  than  the  discarded  coordinate 
direction.  After  the  next  cycle  of  minimizations,  generate  a  new  pattern 
direction  and  again  replace  one  of  the  coordinate  directions.  This  process 
is  illustrated  in  Figure  2. 


Figure  1.  Logic  diagram  for 
interior  penalty  function 
technique 
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20.  Figure  3  Is  a  logic  diagram  for  the  unconstrained  minimization 

algorithm.  The  pattern  move  is  constructed  in  block  A,  then  used  for  a 

minimization  step  (blocks  B  and  C),  and  then  stored  in  (block  D)  as  all 

of  the  directions  are  up-numbered  and  S,  is  discarded.  The  direction  S 

1  n 

will  then  be  used  for  a  minimizing  step  just  before  the  construction  of  the 
next  pattern  direction.  Consequently,  in  the  second  cycle,  both  X  and  Y 
in  block  A  are  points  that  are  minima  along  ,  the  last  pattern  direction. 
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Figure  3.  Logic  diagram  for  minimization 
of  <}i(X) 


This  sequence  will  Impart  special  properties  to  X  -  Y  that  are  the 

source  of  the  rapid  convergence  of  the  method. 

21.  Figure  3  shows  a  block  requiring  a  one-dlmenslonal  minimization 
of  a*  of  the  function  <KX  +  o*S  )  .  The  one-dlmenslonal  minimization 
uses  a  four-point  cubic  Interpolation.  It  finds  the  minimum  along  the 
direction  ,  where  ^  Is  the  coordinate  of  the  previous  minimum.  By 
trial  and  error  It  finds  three  points  with  the  middle  one  less  than  the 
other  two.  It  makes  a  quadratic  Interpolation  and,  then,  a  cubic  Interpola¬ 
tion.  If  the  actual  function  evaluated  at  the  new  Interpolated  point  Is  not 
sufficiently  close  to  that  of  the  preceding  point,  or  If  It  Is  not  sufficiently 
close  to  the  Interpolated  function,  then  another  cubic  Interpolation  Is 
made.  The  logic  for  this  algorithm  Is  shown  In  Figure  4. 
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PART  IV:  APPROXIMATE  COMPUTATION  OF  DOOR  REACTION 


22.  It  should  be  emphasized  that  this  program  is  intended  to  assist 
in  rapid  approximate  design  and  not  detailed  analysis.  The  basic  procedures 
in  TM  5-1300,  TN-1494  and  used  herein  have  been  found  to  be  sufficiently 
accurate  for  simple  geometries  of  beams  and  slabs  without  openings.  Figures 
5  and  6  compare  deflections  for  a  plate  fixed  on  four  sides  and  for  a  beam; 
the  approximate  solutions  and  the  finite  element  solutions  agree  within 
about  10  percent.  However,  Figure  7  shows  that  the  static  shear  procedures 
suggested  in  TM  5-1300  are  substantially  below  dynamic  shears;  this  is  a 
limitation  of  the  approximate  procedures  and  is  under  current  investigation. 

23.  A  steel  door  attached  to  a  concrete  wall  was  examined  using  a 
finite  element  technique.  Figure  8  shows  the  slab  and  door;  Figure  9  shows 
the  deflection  of  the  door  by  the  approximate  procedure  developed  herein  and 
the  finite  element  procedure.  There  is  some  disagreement  in  deflection, 
especially  when  one  considers  the  deflecting  top  support.  It  should  be 
particularly  noted  that  the  deflecting  support  condition  for  actual  doors  on 
slabs  (modeled  correctly  by  finite  element  and  assumed  rigid  by  approximate 
solution)  absorbs  significant  amounts  of  energy  by  rigid-body /door  motion. 
Thus,  the  resulting  center  door  deflection  is  reduced.  The  resulting  dynamic 
shear  around  the  door  (transferred  to  the  wall)  is  reduced  from  what  would 

be  computed  for  a  nondeflecting  plate  using  approximate  dynamic  plate  theory 
(Figure  10).  The  alternatives  are  to  use  finite  element  analysis  procedures 
or  to  modify  dynamic  plate  theory.  Finite  element  analysis  is  certainly  the 
better  approach;  however,  it  is  basically  an  analysis  technique  and  is  more 
difficult  and  expensive  to  use  than  the  simpler  approximate  procedure.  It 
is  suggested  that  the  shear  calculated  from  approximate  plate  theory  be 
adjusted  by  a  constant  for  use  as  a  door  reaction  required  for  input  to  wall 
design  (Ferritto  1977). 

24.  The  maximum  reaction  (REA)  occurs  at  the  moment  the  slab  first 
reaches  yield.  At  this  point  the  combination  of  load  and  resistance  is 
maximum.  Table  1  gives  maximum  dynamic  reaction  for  a  simply  supported 
plate.  For  the  case  of  one  side  free  and  three  sides  simply-supported,  the 
b-dlmenslon  doubled  may  be  used.  The  values  of  pressure  F  and  resistance 
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Figure  5.  Displacement  history  of  4-  by  4-ft 
(1.2-  by  1.2-m)  plate 

R  should  correspond  to  time  of  yielding.  The  reaction  values  should  be 
adjusted  for  support  deflection.  The  value  of  1.0  Is  suggested  for  non¬ 
deflecting  supports  and  0.5  for  full  deflecting  supports  as  approximate 
factors.  Once  design  has  been  finalized  it  is  suggested  that  results  be 
analyzed  using  a  finite  element  analysis. 


DEFLECTION,  IN. 


r 


0  10  20  30 

TIME,  MS 

Figure  6,  Deflection  of  center,  10-ft  (3-m)  beam 


Figure  7.  Shear  in  plate 


CONCRETE  WALL  =  7.125  IN 
FREE  EDGE 


12  IN. 


STEEL  DOOR 


t  =  0.3125  IN. 


Table  1 

Four  Sides,  Uniform  Load* 


Strain 

Range 

Elastic 


Plastic 


SIMPLE 

SUPPORT 

s 

SHORT  sioe 

b 

LONG  Sloe 


a/b 


_ Dynamic  Reactions** 

VA/b 


1.0 

o 

b 

+ 

0.18R 

0.07P 

+ 

0.18R 

0.9 

0.06P 

+ 

0.16R 

0.08P 

+ 

0.20R 

0.8 

0.06P 

+ 

0.14R 

0.08P 

+ 

0.22R 

0.7 

0.05P 

+ 

0.13R 

0.08P 

+ 

0.24R 

0.6 

0.04P 

+ 

O.llR 

0.09P 

+ 

0.26R 

0.5 

0.04P 

+ 

0.09R 

0.09P 

+ 

0.28R 

1.0 

0.09P 

+ 

0.16R 

m 

0.09P 

+ 

0.16R 

m 

0.9 

0.08P 

+ 

0.15R 

m 

0.09P 

+ 

0.18R 

m 

0.8 

0.07P 

+ 

0.13R 

m 

O.lOP 

+ 

0.20R 

m 

0.7 

0.06P 

+ 

0.12R 

m 

O.lOP 

+ 

0.22R 

m 

0.6 

0.05P 

+ 

O.IOR 

m 

O.lOP 

+ 

0.25R 

m 

0.5 

0.04P 

+ 

0.08R 

m 

O.llP 

+ 

0.27R 

m 

*  Based  on  Information  from  Norris  (1959). 

*  P  *  pressure  at  time  of  yield,  psl 

R  °  elastic  resistance,  psl 

R  >  plastic  resistance,  psl 
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PART  V:  THE  COMPUTER  PROGRAM 


25. 


The  CSDOOR  program  is  composed  of  four  areas: 
a.  Blast  load  determination. 

Structural  analysis  parameters. 

£.  Dynamic  response, 
d.  Optimization. 


Subroutines 


Blast  load  determination 

26.  The  blast-load  determination  is  accomplished  by  subroutines  BLA, 
PIC,  SGRID,  HBA,  RATIO,  GRID,  GAS  INTERP,  EQUIV,  HEDATA,  ARDC,  SHOCK,  and 
TNT.  The  subroutines  read  the  explosive  weight  and  type  and  cell  geometry. 
The  equivalent  spherical  weight  of  TNT  and  the  equivalent  pressure  loading 
is  computed  based  on  the  geometry  of  the  wall  and  charge  location.  Both  the 
shock  pressure  and  its  duration  and  the  gas  pressure  and  its  duration  are 
calculated.  Using  the  duration  and  pressure  data  for  both  shock  and  gas, 
the  program  computes  an  equivalent  triangular  pressure  loading  for  each  part 
and  adds  both  to  produce  the  resultant  shown  in  Figure  11.  The  total  impulse 
is  then  determined. 

Structural  analysis 

27.  The  structural  analysis  is  accomplished  by  subroutines  SSTIFF, 
DOORl,  D00R2,  D00R3,  DOORA,  and  DOORS.  These  routines  compute  the  stiffness, 
resistance,  and  equivalent  mass  of  the  plate  using  input  material  properties 
as  in  TM  5-1300.  Both  flexure  and  shear  are  considered.  Openings  in  plates 
are  allowed. 

Dynamic  response 

28.  The  dynamic  response  calculation  is  accomplished  in  subroutine 
RESP.  The  program  determines  the  response  of  the  plate  modeled  as  an 
equivalent  dynamic  single-degree-of-freedom  system  with  bilinear  stiffness 
and  the  pressure  loading  shown  in  Figure  11.  The  solution  technique  is 
based  on  a  Newmark  iteration  method. 

Optimization 

29.  The  optimization  of  an  initial  design  is  accomplished  in  sub¬ 
routines  OPT,  MINIMZ,  PMINZ,  DMINZ,  GETE,  SUMRY,  TLEFT,  and  GCOMP.  The 
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Figure  11.  Equivalent  pressure  loading 


methodology  used  is  that  of  a  penalty  function  with  individual  minimization 
sequences  being  accomplished  by  the  Powell  method. 

Program  Input 


Data  input  guide 

30.  The  following  sections  describe  the  data  input  phase  of  CSDOOR 
and  the  various  options  available.  A  data  input  guide  was  prepared  to  aid 
the  user  in  data  preparation.  This  guide,  with  appropriate  entries,  is 
presented  in  Appendix  A  with  each  example  problem.  Also,  a  blank  copy  of 
the  guide  is  presented  at  the  back  of  this  report.  Illustrative  results  are 
presented  for  the  following  example  problems: 

a.  Analyze  steel  door. 

b.  Analyze  steel  plate. 

£.  Analyze  steel  wall. 

Data  groups 

31.  Defining  a  problem  Involves  specification  of  7  basic  data  groups 
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composed  of  about  48  variables.  The  program  can  be  run  by  making  use  of  an 
existing  data  file  having  sequence  numbers  at  the  start  of  each  line.  As  an 
alternative  mode  of  Input,  an  Interactive  phase  Is  also  provided  which 
assists  the  user  In  defining  data  for  a  particular  problem.  All  data  are 
entered  In  free-fleld  format  with  commas  or  blanks  used  to  separate  the 
successive  numbers.  All  values  can  be  Input  with  or  without  decimal  points 
(for  Instance,  FLAGl  =  1  can  be  Input  either  as  1.  or  as  1).  If  the  user  so 
desires,  data  Input  Interactively  can  be  saved  Into  a  permanent  file  with 
line  numbers.  The  output  from  a  problem  can  be  written  to  the  terminal  or 
into  a  permanent  file  to  be  either  scanned  with  an  editor  or  sent  to  a  line 
printer. 

32.  The  user  should  be  aware  that  data  saved  in  a  file  may  not 
coincide  exactly  with  the  values  input  interactively.  The  data  are  written 
to  a  file  using  field  widths  adequate  for  practical  situations.  For 
Instance,  most  variables  are  written  using  two  digits  past  the  decimal 
point.  In  the  event  that  greater  accuracy  is  needed  in  the  recorded  data, 
the  data  file  can  be  edited  accordingly. 

33.  The  different  data  groups  with  names  of  the  variables  for  each 
one  as  used  in  the  program  are  as  follows: 

£.  Data  Group  1 — Heading  (HDG) ; 

HDG  -  Alphanumeric  heading  for  problem  identifica¬ 

tion;  68  characters  maximum 

Data  Group  2 — Program  Control  (FLAGl,  FLAG2,  FLAG3, 

FLAG4,  FLAGS.  PC,  COST): 

FLAGl*  -  Set  =  1  for  optimization;  otherwise  =  0 

FLAG2  -  Set  =  0  to  calculate  gas  pressure;  set  =  1  to 

input  gas  pressure 

FLAG3  -  Set  =  0  for  plate  thickness  (TS) ;  set  =  1  for 
section  modulus  (Z)  and  moment  of  inertia  (I) 

FLAG4  -  Set  =  1  for  impulse  grid;  otherwise  =  0 

FLAGS  -  Set  =  1  for  door/window  reaction  present; 

otherwise  =  0 

PC  -  Set  =  0  for  standard  printout 

=  1  for  print  response  time-history 


*  Optimization  cannot  be  used  if  a  composite  door  is  used  (FLAG3  *  1)  or 
if  a  door  is  present  (FLAGS  >  1). 


COST 


-  Cost  of  steel  per  dollars  per  lb*  (default 
•  0.60) 

£.  Data  Group  3 — Load  Parameters  (WLB,  ANUM,  RLOD,  CASE, 

APAMB.  TAMB,  ALTKFT) ; 

WLB  -  Weight  of  actual  explosive  Including  safety 

factor,  lb 

ANUM  -  Explosive  number  used  to  compute  explosive 

equivalence  (see  Table  2  for  list  of  explosives) 

RLOD  -  Explosive  length  to  diameter  ratio  (default  =  1) 

CASE  -  Projectile  case  weight  to  explosive  weight 
ratio  (use  0  for  conservative  analysis) 

APAMB  -  Ambient  air  pressure,  psia  (default  =  14.69) 

TAMB  -  Ambient  temperature,  *C  (default  =  20“C) 

ALTKFT  -  Altitude,  10^  ft  (when  APAMB  and  TAMB  not 
specified) 

Data  Group  4 — Geometry: 

(1)  When  gas  pressure  is  calculated  (FLAG2  =  0)  input 
(RA,  H,  EL,  HLIT,  ELLIT,  AV.  AC,  ICODE(i).  where 
1  =  1 ,  2 ,  3 ,  and  4 ) ; 


RA  -  Distance  from  charge  to  wall,  ft 

H  -  Wall  height,  ft 

EL  -  Wall  length,  ft 

HLIT  -  Height  of  charge,  ft 

ELLIT  -  Distance  of  charge  to  left  boundary,  ft 

3 

AV  -  Cell  volume  for  gas  pressure,  ft 

2 

AC**  -  Cell  vent  area  for  gas  pressure,  ft 

ICODE(l)  -  Set  =  1  for  floor  reflection;  otherwise  set 
=  0 


IC0DE(2)  -  Set 
=  0 


1  for  roof  reflection;  otherwise  set 


IC0DE(3)  -  Set  =  1  for  left  wall  reflection;  otherwise 
set  =  0 

IC0DE(4)  -  Set  =  1  for  right  wall  reflection;  otherwise 
set  =  0 


A  table  of  factors  for  converting  U.  S.  customary  (NON-SI)  units  of  mea¬ 
surement  to  metric  (SI)  units  is  presented  on  page  5. 

CSDOOR  will  not  solve  for  gas  pressure  if  vent  area  ■  0. 


Table  2 

List  of  Explosives 


Explosive 

Ntunber 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 
29 


Explosive  Name  and  Composition 
TNT 
TNETB 

EXPLOSIVE  D 

PENTOLITE  (PETN/TNT  50/50) 
PICRATOL  (EXPLOSIVE  D/TNT  52/48) 
CYCLOTOL  (RDX/TNT  70/30) 

COMP  B  (RDX/TNT/WAX  59.4/39.6/1.0) 
RDX/WAX  (98/2) 

COMP  A-3  (RDX/WAX  91/9) 
TNETB/AL  (90/10) 

TNETB/AL  (78/22) 

TNETB/AL  (72/28) 

TNETB/AL  (65/34) 

TRITONAL  (TNT/AL80/70) 
RDX/AL/WAX  (88/10/2) 
RDX/AL/WAX  (89/20/2) 
RDX/AL/WAX  (74/21/5) 
RDX/AL/WAX  (74/22/4) 
RDX/AL/WAX  (62/33/5) 

TORPEX  II  (RDX/TNT/AL  42/40/18) 
H6  (RDX/TNT/AL/WAX  45/29/21/5) 
HBX-1  (RDX/TNT/AL/WAX  40/38/16/5) 
HBX-3  (RDX/TNT/AL/WAX  31/29/35/5) 
TNETB/RDX/AL  (39/26/35) 
ALUMINUM 
WAX 
RDX 
PETN 
TETRYL 
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I 


e. 


f. 


(2)  When  gas  pressure  Is  Input  (FLAG2  »  1)  input  (TOTIM, 
H,  EL,  FPRES,  TO,  PG.  TG,  ICODE(i).  where  i  =  1,  2, 
3.  and  4); 

TOTIM  -  Total  impulse,  psi-msec 

H  -  Wall  height,  ft 

EL  -  Wall  length,  ft 

FPRES  -  Peak  pressure,  psi 

TO  -  Duration  of  peak  pressure,  msec 

PG  -  Gas  pressure,  psi 

TG  -  Gas  pressure  duration,  msec 

ICODE(l)  -  Set  =  1  for  floor  reflection;  otherwise 
set  =  0 

IC0DE(2)  -  Set  =  1  for  roof  reflection;  otherwise 
set  =  0 

IC0DE(3)  -  Set  =  1  for  left  wall  reflection;  otherwise 
set  =  0 

IC0DE(4)  -  Set  =  1  for  right  wall  reflection;  otherwise 
set  =  0 


Data  Group  5 — Strength  Parameters  (FDY,  TS,  SN,  DH,  DEL, 

U.  E); 

FDY  -  Steel  dynamic  design  strength,  psi 

TS*  -  Door  thickness,  in. 

SN  -  Support  code  (see  Figure  12) : 

=  1,  bottom  fixed 
=  2,  bottom  and  1  side  fixed 
=  3,  bottom  and  2  sides  fixed 
=  4,  4  sides  fixed 

=  5,  beam  simple  supports  top  and  bottom 
=  6,  beam  fixed  top  and  bottom 
=  7,  beam,  simple  support  top,  fixed  bottom 
=  13,  three  sides  simple  supports,  bottom  free 
=  14 ,  four  sides  simple  support 

DH*  -  Door  height,  ft 

DEL*  -  Door  width,  ft 

U  -  Allowable  ductility  limit  for  optimization 

E  -  Modulus  of  elasticity,  psi 

Data  Group  6 — Door  Properties  (ZHOR,  ZVERT,  lAVG,  WDR) ; 

ZHOR  -  Plastic  horizontal  section  modulus/in. , 

in.^/in. 


*  Door  and  wall  are  s3mon}nnous  when  door  height  and  width  equal  wall 
height  and  width. 
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3 

ZVERT  -  Plastic  vertical  section  modulus/in.,  in.  /in. 

4 

lAVG  -  Average  moment  of  inertia/in.,  in.  /in. 

WDR  -  Total  door  weight,  lb 

£.  Data  Group  7 — Opening  (Panel)  Parameters*  (see  Fig 

ure  13)  when  FLAGS  =  1  input  (H2.  WT,  B,  REA,  RDl,  Hi): 

H2  -  Opening  height,  ft 

WT  -  Opening  width,  ft 

B  -  Distance  from  left  side  to  opening,  ft 

REA  -  Opening  reaction,  Ib/in.  (3  sides  supported) 

RDl  -  Resistance  for  calculating  opening  reaction, 

lb  (3  sides  supported) 

HI  -  Distance  to  floor,  ft  (for  window  only) 

34.  The  explosive  number  (Data  Group  3)  refers  to  the  list  of  explo¬ 
sives  in  Table  2.  This  number  is  used  to  compute  explosive  equivalence. 

The  length/diameter  ratio  for  an  explosive  sphere  is  0.0,  which  gives  a 
shape  factor  of  1.0.  For  an  uncased  explosive  the  case  explosive  weight 
ratio  is  0.0.  For  sea  level  calculations,  the  ambient  air  pressure  > 

temperature  altitude  can  be  left  blank  and  will  default  to  14.69 

psi  and  20®C.  If  FLAG2  in  Data  Group  2  is  set  to  1,  the  impulse,  duration, 
and  pressure  will  be  read  on  Date  Group  4.  If  FLAG2  is  set  to  zero,  the 
charge  to  wall  distance,  charge  height,  and  distance  from  the  left  side  will 
be  read.  If  SN  is  zero,  the  program  will  sum  the  number  of  reflecting 
sidewall  surfaces  specified  in  Data  Group  4.  The  separate  use  of  SN  is 
helpful  when  a  frangible  wall  is  present,  which  creates  a  shock  reflection 
but  does  not  provide  any  support. 

35.  The  SN  (see  Figure  12)  conditions  1  through  4  are  intended  to  be 
used  to  represent  steel  cell  walls  and  roofs;  SN  conditions  5  through  7  are 
steel  plates  spanning  in  one  direction.  The  SN  conditions  13  and  14  are 
specifically  intended  to  represent  typical  steel  plate  doors  and  pass¬ 
through  windows. 


*  When  an  opening  (panel)  is  present,  the  program  analyzes  the  wall  as  a 
door.  The  panel  reaction  must  be  provided  either  as  a  resistance  (Ib/ln.) 
along  the  edges  of  the  panel  or  as  a  resistance  (psi)  over  the  entire 
panel.  This  case  seldom  occurs  in  practice. 
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•NOTE:  OPENING  MUST  BE  IN  CENTER  OF  WALL. 


WALL  FOUR  SIDES  SUPPORTED  WITH  OPENING 


Figure  13.  Wall  geometry  with  opening 


Example  Problems 

36.  Three  example  problems  are  presented  in  Appendix  A.  In  examples  1 
and  2,  data  were  entered  from  a  data  file.  In  example  3,  data  were  entered 
Interactively. 


i 
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APPENDIX  A:  EXAMPLE  PROBLEM  1,  ANALYZE  STEEL  DOOR 


Given 

Geometry  as  shown 

Charge  Wt.  =  10.0  lb  Comp  B 

(uncased) 

Floor,  roof,  left  wall  and 
right  wall  deflection 

FDY  =  48000  psi 

Three  sides  simply  supported 
bottom  free 

Allowable  ductility  y  =  3 


Required 


Analyze  steel  door 


Assume 

3 

Plastic  section  modulus  horiz.  =  1.77  in.  /In. 

3 

Plastic  section  modulus  vert.  =  1.51  in.  /in. 
Average  moment  of  inertia  =  3.495  in.^/in. 


Door  weight  =  920  lbs 


A1 


STEEL  DOOR 


ELEVATION 


^  STEEL  DOOR 
Z! — yLLJiJU  I 


10# COMP. 


I  3x5  I 

ROOF  VENT  I 

J  ' 


FILENAME  DOOREX1 


optlodtatLon  cannot  b«  um4  If  fLAC3 
Door  and  wail  ara  aynonyaeus  when  OM 
FLM;)  •  1  can  only  ba  uaad  whan  SN  * 


1000  EXAMPLE  PROBLEM  ONE 

1010  0010000 

1020  10  7  1  0  0  0  0 

1030  3  10  12  3  5  1200  15  1  1  1  1 

1040  48000  0  13  6  4  3  0 

1050  1.77  1.51  3.495  920 

0X0057 

THIS  PROGRAM  SHOULD  BE  USED  ONLY  BY  ENGINEERS 
WHO  ARE  EXPERIENCED  IN  BLAST  DESIGN  AND  ARE 
THOROUGHLY  FAMILIAR  WITH  METHODS  OF  ANALYSIS 
DESCRIBED  IN  TM  5-1300  'STRUCTURES  TO  RESIST 
THE  EFFECTS  OF  ACCIDENTAL  EXPLOSIONS.'  CONNECTIONS 
AND  DETAILS  MUST  BE  CAREFULLY  DESIGNED  TO  ACHIEVE 
THE  DEGREE  OF  FIXITY  THAT  IS  ASSUMED  IN  THE  PROGRAM 

INPUT  NAME  OF  DATA  FILE  IN  7  CHARACTERS  OR  LESS 

HIT  A  CARRIAGE  RETURN  IF  DATA  TO  COME  FROM  TERMINAL. 

I>D00REX1 

INPUT  NAME  OF  FILE  FOR  OUTPUT  TO  BE  WRITTEN  TO. 

HIT  A  CARRIAGE  RETURN  IF  OUTPUT  TO  TERMINAL 

I> 


EXAMPLl'.  PROBLEM  ONE 
COMP  B  <RDX/TNT/UAX»59. 4/39. 6/1.0) 

EXPLOSIVE  PROPERTIES . CHARGE  UEIGHT(LB)  =  10.00 

NUMBER  EQUT  EFORM  EXPLOSIVE  COMPOSITION  BY  WEIGHT 
KCAL/G  C  H  N  0  AL 

7  1.100  .004330  .252  .026  .298  .424  0.000 

PAMB(PSIA)=  14.69  TAMB<C)=  20.00 

SHOCK  WAVE  CALCULATION 

INPUT  PARAMETERS  CHARGE  WEIGHT  ADJUSTMENTS 

CHARGE  WEIGHT (LB)  =  10.00  ADJUSTED  WT(LB  TNT)  =  1 

EXPLOSIVE  NUMBER  =7  HE  ENERGY  FACTOR  =  1 

L/D  RATIO  =  1.000  CHARGE  SHAPE  FACTOR  =  1 

CASE/CHARGE  WT  RATIO  *  0.  CASE  HEIGHT  FACTOR  =  1 


.00 

100 

000 

000 


I 


r 


CHAMBER  PRESSURE (PSIA>b  14.69 
CHAMBER  TEHP(C)  s  20.00 
ALTITUDE  (KFT)  »  0. 

DESIRED  DISTANCE  (FT)  s  3.000 
(CM)  =  91.44 


PRESSURE  SCALE  FACTOR^  1.000 
DISTANCE  SCALE  FACTOR-  .4496 
TINE  SCALE  FACTOR  *  .4535 
NORMAL  REFL'  FACTOR  >  7.526 


TINE  AFTER 

TINE  AFTER 

INCIDENT 

NORM  REFL 

EXPLOSION 

SHOCK  ARR 

OVERPRESS 

OVERPRESS 

(MSEC) 

(MSEC) 

(PSI) 

(PSI) 

.2615 

0. 

497.2 

3742. 

.3698 

.1083 

156.8 

1180. 

.4240 

.1625 

98.81 

743.6 

.4782 

.2167 

64.16 

482.9 

.5324 

.2708 

41.97 

315.9 

.5865 

.3250 

27.05 

203.6 

.6407 

.3792 

16.66 

125.4 

.6949 

.4334 

9.259 

69.68 

.7490 

.4875 

3.905 

29.39 

.8032 

.5417 

0. 

0. 

IMPULSE  (PSI.NSEO  — 

INCIDENT  -  53.86 

REFLECTED=  405.3 

. CAUTION— CNTACT  SURFACE  HAS  ARRIVED. 

DATA  ARE  CRUDE  BEYOND  T(MSEC>  AFTER  SHOCK  ARRIVAL^^ 


7.5617E-02 


DISTANCE  OF  CHARGE  FROM  BLAST  WALL 

FT. 

3.00 

CHARGE  WEIGHT 

LBS. 

11.00 

BLAST  WALL  HEIGHT 

FT. 

10.00 

BLAST  WALL  LENGTH 

FT. 

12.00 

HEIGHT  OF  CHARGE  ABOVE  GROUND 

FT. 

3.00 

MIN.  DIST.  BETWEEN  CHARGE  *  ADJ.  HALL 
REFLECTION  CODE 

FT. 

5.00 

1111 

TOTAL  IMPULSE  440.28  PSI-MS 

VENT  AREA  15.00  FT2  CELL  VOLUME  1200.00  FT3 


OAS  PRESSURES  CALCULATION 


PEAK  GAS  PRESSURE 
GAS  DURATION 
GAS  IMPULSE 
TOTAL  IMPULSE 


68.85  PSI 
91.84  MSEC 
3161.46  PSI-MSEC 
3299.79  PSI-MSEC 


A5 


DURATION  OF  LOAD 


S.721S6  NSEC 


FICTITIOUS  PEAK  PRESSURE  153.90  PSI 


EFFECTIVE  IHPULSE  3299.79  PSI-HSEC 


FS  DYNAMIC 

48000.00 

PSI 

PLATE  THICKNESS 

0.00 

IN 

SUPPORT  CODE 

13.00 

DOOR  HEIGHT 

6.00 

FT 

DOOR  LENGTH 

4.00 

FT 

PLASTICITY  (MU) 

3.00 

HORIZONTAL  Z 

1.77 

IN3/IN 

VERTICAL  Z 

1.51 

IN3/IN 

I  AVERAGE 

3.50 

IN4/IN 

DOOR  WEIGHT 

920.00 

LBS 

HEIGHT 

72.00  IN 

LENGTH 

48.00  IN 

POSITIVE 

NEGATIVE 

POSITIVE 

NEGATIVE 

VERTICAL  MOMENT 
VERTICAL  MOMENT 
HORIZONTAL  MOMENT 
HORIZONTAL  MOMENT 

72480.00 

72480.00 

84960.00 

84960.00 

IN-LBS/IN 

IN-LBS/IN 

IN-LBS/IN 

IN-LBS/IN 

WIDTH 

WIDTH 

WIDTH 

WIDTH 

THREE  SIDES  SIHPLY  SUPPORTED 


X 

24.0000 

IN 

Y 

32.7453 

IN 

RU 

417.4331 

PSI 

W1 

423.3625 

PSI 

W2 

405.5744 

PSI 

XE 

.2909 

IN 

K 

1435.13 

LB/IN/IN2 

HASS 

469.15 

LB-MSEC2/IN/IN2 

ALLOWABLE  HAX  DEFLECTION 


.8726  IN 


HASS 

LOAD 

DURATION 

RESISTANCE 


469.150  LB-HSEC2/IN/IN2 
153.902  PSI 
5.722  MSEC 
417.433  PSI 


STIFFNESS 

GAS  PRESSURE 
DURATION 

1435.131 

68.85 

91.84 

LB/IN/IN2 

PSI 

MSEC 

MEMBRANE  YIELD  DEFLECTION 

1.518582 

IN 

ELASTIC  DEFLECTION  LIMIT 

.290868 

IN 

MAXIMUM  DEFLECTION 

.181843 

IN 

NATURAL  PERIOD 

3.592439 

MSEC 

TIME  TO  MAXIMUM 

DEFLECTION 

1.678216 

MSEC 

DURATION/NATURAL 

PERIOD 

25.565389 

LOAD/RESISTANCE 

.368687 

1.200  CP  SECONDS  EXECUTION  TINE. 
C> 


APPENDIX  B:  EXAMPLE  PROBLEM  2,  ANALYZE  STEEL  PLATE 

*  ► 


Given 

Geometry  as  shown 
Charge  wt  =  10  lbs  T.N.T. 

Explosive  length  to  diameter  ratio  2.5 
Case  to  explosive  ratio  1.2 
Floor,  left  wall  and  right  wall  reflection 
Four  sides  simply  supported 
FDY  =  48000  psi 
Allowable  ductility  p  =  5 
Required 

Analyze  steel  plate 

Assume 

Plate  thickness  0.75  in. 


*■ 

1 


k 


{ 


I 

t 


c 


B1 


DOOREX2 


T 


1000  EXAMPLE  PROBLEM  TWO 

1010  0000000 

1020  10  1  2.5  1.2  14.69  20  0 

1030  4  10  12  35001011 

1040  48000  0.75  14  1.33  4  5  0 

0X0057 

THIS  PROGRAM  SHOULD  BE  USED  ONLY  BY  ENGINEERS 
WHO  ARE  EXPERIENCED  IN  BLAST  DESIGN  AND  ARE 
THOROUGHLY  FAMILIAR  WITH  METHODS  OF  ANALYSIS 
DESCRIBED  IN  TH  5-1300  'STRUCTURES  TO  RESIST 
THE  EFFECTS  OF  ACCIDENTAL  EXPLOSIONS.'  CONNECTIONS 
AND  DETAILS  MUST  BE  CAREFULLY  DESIGNED  TO  ACHIEVE 
THE  DEGREE  OF  FIXITY  THAT  IS  ASSUMED  IN  THE  PROGRAM 

«««««««t««««««««««««CAUTION«««««««««««*«»*«««** 

INPUT  NAME  OF  DATA  FILE  IN  7  CHARACTERS  OR  LESS 

HIT  A  CARRIAGE  RETURN  IF  DATA  TO  COME  FROM  TERMINAL. 

I>D00REX2 

INPUT  NAME  OF  FILE  FOR  OUTPUT  TO  BE  WRITTEN  TO. 

HIT  A  CARRIAGE  RETURN  IF  OUTPUT  TO  TERMINAL 
I> 


EXAMPLE  PROBLEM  TWO 

TNT 

EXPLOSIVE  PROPERTIES . CHARGE  WEIGHKLB)  =  10.00 

NUMBER  EQHT  EFORM  EXPLOSIVE  COMPOSITION  BY  HEIGHT 
KCAL/6  C  H  N  0  AL 
1  1.000  -.078400  .370  .022  .185  .423  0.000 


PAHB(PSIA)*  14.69  TAHB(C)=  20.00 

. CHARGE  SHARE  CORRECTION  IS  CRUDE.  PSI  EXCEEDS  RANGE  OF  EXPERIMENTAL  DATA. 

. CASE  HEIGHT  CORRECTION  IS  CRUDE.  PSI  EXCEEDS  RANGE  OF  EXPERIMENTAL  DATA. 


INPUT  PARAMETERS 

SHOCK  HAVE 

CALCULATION 

CHARGE  HEIGHT  ADJUSTMENTS 

CHARGE  HEIGHT (LB) 
EXPLOSIVE  NUMBER 

L/D  RATIO 

»  10.00 
=  1 
>  2.500 

ADJUSTED  UT(LB  TNT)  = 

HE  ENERGY  FACTOR 

CHARGE  SHAPE  FACTOR  <= 

21.41 

1.000 

3.012 

CASE/CHARGE  UT  RATIO 

m 

1.200 

CASE  HEIGHT  FACTOR 

.7109 

CHAMBER  PRESSURE <PS1 A) 

m 

14.69 

PRESSURE  SCALE  FACTOR* 

1.000 

CHAMBER  TEMP(C) 

s 

20.00 

DISTANCE  SCALE  FACTOR* 

.3601 

ALTITUDE  (KFT) 

a 

0. 

TIME  SCALE  FACTOR  - 

.3632 

DESIRED  DISTANCE  (FT) 

s 

4.000 

NORMAL  REFL*  FACTOR  » 

7.307 

(CM)  * 

121.9 

TIME  AFTER 

TIME  AFTER 

INCIDENT 

NORN  REFL 

EXPLOSION 

SHOCK  ARR 

OVERPRESS 

OVERPRESS 

(MSEC) 

(MSEC) 

(PSD 

(PSD 

.3687 

0. 

437.8 

3199. 

.5129 

.1442 

138.1 

1009. 

.5850 

.2163 

87.00 

635.7 

.6571 

.2884 

56.49 

412.8 

.7292 

.3605 

36.96 

270.0 

.8013 

.4326 

23.82 

174.0 

.8733 

.5047 

14.67 

107.2 

.9454 

.5767 

8.152 

59.57 

1.018 

.6488 

3.438 

25.12 

1.090 

.7209 

0. 

0. 

IHPULSE  (PSI.HSEO  — 

INCIDENT  »  63.11 

REFLECTED-  461.1 

. CAUT10N--CNTACT  SURFACE  HAS  ARRIVED. 

DATA  ARE  CRUDE  BEYOND  T(HSEC>  AFTER  SHOCK  ARRIVAL-  .1163 


DISTANCE  OF  CHARGE  FROM  BLAST  HALL 

FT. 

4.00 

CHARGE  HEIGHT 

LBS. 

21.41 

BLAST  HALL  HEIGHT 

FT. 

10.00 

BLAST  HALL  LENGTH 

FT. 

12.00 

HEIGHT  OF  CHARGE  ABOVE  GROUND 

FT. 

3.00 

MIN.  BIST.  BETUEEN  CHARGE  F  ADJ. 

HALL 

FT. 

5.00 

REFLECTION  CODE 

1  0 

1  1 

TOTAL  IHPULSE 

542.80 

PSI-MS 

DURATION  OF  LOAD 

5.04033 

MSEC 

FICTITIOUS  PEAK  PRESSURE 

215.38 

PSI 

EFFECTIVE  IMPULSE 

542.80 

PSI-MSEC 

F8  DYNAMIC  48000.00  P8I 


PLATE  THICKNESS 

.75 

IN 

SUPPORT  CODE 

14.00 

DOOR  HEIGHT 

1.33 

FT 

DOOR  LENGTH 

4.00 

FT 

PLASTICITY  (HU) 

5.00 

HEIGHT  15.96  IH  LEHGTH 


46.00  IN 


POSITIVE  VERTICAL  HQHENT 
NEGATIVE  VERTICAL  HOHENT 
POSITIVE  HORIZONTAL  HOHENT 
NEGATIVE  HORIZONTAL  HOHENT 


6750.00  IN-LSS/IN  WIDTH 
6750.00  IN-LBS/IN  WIDTH 
6750.00  IN-LBS/IN  WIDTH 
6750.00  IN-LBS/IN  WIDTH 


FOUR  SIDES  SIHPLY  SUPPORTED 


X  11.42  IN 

T  7.9B  IN 

RU  310.50  PSI 


XE  .2165  IN 

K  1434.05  LB/1N/IN2 

HASS  379.42  LB-HSEC2/1N/1N2 


ALLOWABLE  HAX  DEFLECTION 


1.0826  IN 


HASS 

LOAD 

DURATION 
RESISTANCE 
STIFFNESS 
GAS  PRESSURE 
DURATION 


379.423 

215.382 

5.040 

310.501 

1434.052 

0.00 

0.00 


LB-HSEC2/IN/IN2 

PSI 

HSEC 

PSI 

LB/1N/IN2 

PSI 

HSEC 


HEHBRANE 

YIELD  DEFLECTION 

.504928 

IN 

ELASTIC 

DEFLECTION  LIHIT 

.216520 

IN 

HAXIHUH 

DEFLECTION 

.260201 

IN 

NATURAL 

PERIOD 

3.231908 

HSEC 

TIHE  TO  HAXIHUH  DEFLECTION 
DURATION/NATURAL  PERIOD 
LOAD/RESISTANCE 
CALCULATED  DIF  FOR  FDY 
TIHE  TO  YIELD 


1.588217  HSEC 
1.559551 
.693661 
1.500000 

1.13059537  HSEC 


1.434  CP  SECONDS  EXECUTION  TIHE. 
C> 


I 


APPENDIX  C:  EXAMPLE  PROBLEM  3,  ANALYZE  STEEL  WALL 


Given 

Geometry  as  shovm 
Pressure-time  history 
Bottom  and  two  sides  fixed 
Allowable  ductility  y  =  3 
No  reflecting  surfaces 
FY  =  48000  psi 
Required 

Analyze  steel  wall 

Assume 

Wall  properties  as  given 

Note 

Data  input  from  terminal 


Cl 


WALL  ELEVATION 


WALL  PROPERTIES:  Ig  =  9.72  IN.^/IN. 

Zhor  =3.67  IN.^/IN. 
ZvERT  =  ’-^2  IN.=*/IN. 
WDR  =  7368  LB 


C2 


UADIttC  (JOt  OBSCau 


THIS  PROGRAM  SHOULD  BE  USED  ONLY  BY  ENGINEERS 
UHO  ARE  EXPERIENCED  IN  BLAST  DESIGN  AND  ARE 
THOROUGHLY  FAMILIAR  WITH  METHODS  OF  ANALYSIS 
DESCRIBED  IN  TM  5-1300  'STRUCTURES  TO  RESIST 
THE  EFFECTS  OF  ACCIDENTAL  EXPLOSIONS.'  CONNECTIONS 
AND  DETAILS  MUST  BE  CAREFULLY  DESIGNED  TO  ACHIEME 
THE  DEGREE  OF  FIXITY  THAT  IS  ASSUMED  IN  THE  PROGRAM 

INPUT  NAME  OF  DATA  FILE  IN  7  CHARACTERS  OR  LESS 
HIT  A  CARRIAGE  RETURN  IF  DATA  TO  COME  FROM  TERMINAL. 

I> 


ENTER  CONUERSIONAL  MODE  FOR  DATA  INPUT 


INPUT  NAME  OF  FILE  DATA  IS  TO  BE  WRITTEN  TO. 

HIT  A  CARRIAGE  RETURN  IF  YOU  DO  NOT  WANT  THIS  FILE. 

I>D00REX3 

INPUT  A  QUESTION  MARK  (?)  IF  MORE  INFORMATION  IS  NEEDED 
INPUT  HEADING  (HDG): 

I>  EXAMPLE  PROBLEM  THREE 

INPUT  PROGRAM  CONTROL  (FLAGl  fFLAG2»FLAG3»FLAG4>FLAG5rPC>C0ST>  *. 

I>0>lf If0>0t0>0 

INPUT  GEOMETRY  ( TOTIMrHf EL »FPRES»TO»PC»TG» ICODE ( I ) »  WHERE  I=1.2»3»4): 
I>4126f lOf 12f343.2f 4.69?110.3f68.1r0r0>0>0 

INPUT  STRENGTH  PARAMETERS  ( FDY  rTS  r  SNf  DHf  DEL  r  MU  > EC )  *. 

I  48000  »0»3»0»0»3f0 

PANEL  PROPERTIES  < ZHOR> ZVER> lAVG* WDR) : 

I>3.&7.1.42»9.72r73A8 

INPUT  NAME  OF  FILE  FOR  OUTPUT  TO  BE  WRITTEN  TO. 

HIT  A  CARRIAGE  RETURN  IF  OUTPUT  TO  TERMINAL 
I>DOOROT3 

.086  CP  SECONDS  EXECUTION  TIME. 


C>0LD»D00R0T3 

OLIST 


EXAMPLE 

PROBLEM  TH 

REE 

BLAST  HALL 

HEIGHT  10.00 

FT 

BLAST  WALL 

LENGTH  12.00 

FT 

DURATION  OF 

LOAD  4.69000 

MSEC 

FICTITIOUS 

PEAK  PRESSURE 

343.20 

EFFECTIVE  IHPULSE  412<i.C>0 


PSI 

PSI-HSEC 


FS  DYNAMIC 
PLATE  THICKNESS 
SUPPORT  CODE 
DOOR  HEIGHT 
DOOR  LENGTH 
PLASTICITY  <HU) 
HORIZONTAL  Z 
VERTICAL  Z 
1  AVERAGE 
DOOR  HEIGHT 


48000.00 

PSI 

0.00 

IN 

3.00 

0.00 

FT 

0.00 

FT 

3.00 

3.67 

IN3/IN 

1.42 

IN3/IN 

9,72 

IN4/IN 

7368.00 

LBS 

HEIGHT  120.00  IN  LENGTH 


144.00  IN 


POSITIVE  VERTICAL  MOMENT 
NEGATIVE  VERTICAL  MOMENT 
POSITIVE  HORIZONTAL  MOMENT 
NEGATIVE  HORIZONTAL  MOMENT 


48160.00  IN-LBS/IN  WIDTH 
68160.00  IN-LBS/IN  WIDTH 
176160.00  IN-LBS/IN  WIDTH 
176160.00  IN-LBS/IN  WIDTH 


SUPPORT  ON  3  SIDES 


LOCATION  YIELD  LINE  LENGTH 
LOCATION  YIELD  LINE  HEIGHT 
ULTIMATE  LOAD  CAPACITY  RU 
SHEAR  LOAD  AT  HORIZ  SUPPORT 
SHEAR  LOAD  AT  VER  SUPPORT 


72.00  IN 
60.69  IN 
185.0344  PSI 
10869.66  LB/IN  WIDTH 
6738.18  PSI 


C5 


LOAD  HASS  FACTOR 
HASS 


.6548 

724.79  LB-HSEC2/IM/1N2 


FIRST  YIELD  POINT  AT  PT3 
ELASTIC  LIHIT  RE  PSI 

ELASTIC  DEFLECTION  XE 

65.20 

.2967 

IN 

SECOND  YIELD  AT  PT  2 

ELASTO  PLASTIC  LIHIT 

93.35 

PSI 

ELASTO-PLASTIC  DEFLECTION 

.6349 

IN 

ULTIHATE  RESISTANCE 

185.03 

PSI 

PLASTIC  DEFLECTION 

2.1888 

IN 

ULTIHATE  RESISTANCE  RU 

185.03 

PSI 

ELASTIC  DEFLECTION  LIHIT  XE 

1.6455 

IN 

STIFFNESS  KE 

112.45 

LB/ 

ALLOWABLE  HAX  DEFLECTION 

4.9364 

IN 

HASS 

LOAD 

DURATION 
RESISTANCE 
STIFFNESS 
GAS  PRESSURE 
DURATION 


724.795  LB-HSEC2/IN/IN2 
343.200  PSI 
4.690  HSEC 
185.034  PSI 
112.450  LB/1N/IN2 
110.30  PSI 
68,10  HSEC 


HEHBRANE 

YIELD  DEFLECTION 

3.796454 

IN 

ELASTIC 

DEFLECTION  LIHIT 

1.645479 

IN 

HAXIHUH 

DEFLECTION 

3.665975 

IN 

NATURAL 

PERIOD 

15.951698 

HSEC 

TIHE  TO 

HAXIHUH  DEFLECTION 

9.037378 

HSEC 

DURATION/NATURAL  PERIOD 

4.269138 

C6 


r 


LOAD/RESISTANCE 
CALCULATED  DIF  FOR  FDY 
TIHE  TO  YIELD 


1.8S4790 

1.447980 

3.28172933 


C> 


MSEC 


C7 


APPENDIX  D;  SAMPLE  DATE  GUIDE 


i 

I 

I 

1 

c 
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Instruction  Report  K-80-3 
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Technical  Report  K-80-4 
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Instruction  Report  K-81-7 
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Technical  Report  K-81-2 
Instruction  Report  K-82-6 
Instruction  Report  K-82-7 


Title 

List  of  Computer  Programs  for  Computer-Aided  Structural  Engineering 

User’s  Guide:  Computer  Program  with  Interactive  Graphics  for 
Analysis  of  Plane  Frame  Structures  (CFRAME) 

Survey  of  Bridge-Oriented  Design  Software 

Evaluation  of  Computer  Programs  for  the  Design/Analysis  of 
Highway  and  Railway  Bridges 

User's  Guide:  Computer  Program  for  Design/Review  of  Curvi¬ 
linear  Conduits/Culverts  (CURCON) 

A  Three-Dimensional  Finite  Element  Data  Edit  Program 

A  Three-Dimensional  Stability  Analysis/Design  Program  (3DSAD) 
Report  1:  General  Geometry  Module 
Report  3:  General  Analysis  Module  (CGAM) 

Report  4:  Special-Purpose  Modules  for  Dams  (CDAMS) 

Basic  User's  Guide:  Computer  Program  for  Design  and  Analysis 
of  Inverted-T  Retaining  Walls  and  Floodwalls  (TWDA) 

User's  Reference  Manual:  Computer  Program  for  Design  and 
Analysis  of  Inverted-T  Retaining  Walls  and  Floodwalls  (TWDA) 

Documentation  of  Finite  Element  Analyses 
Report  1:  Longview  Outlet  Works  Conduit 
Report  2.  Anchored  Wall  Monolith.  Bay  Springs  Lock 

Basic  Pile  Group  Behavior 

User's  Guide.  Computer  Program  for  Design  and  Analysis  of  Sheet 
Pile  Walls  by  Classical  Methods  (CSHTWAL) 

Report  1:  Computational  Processes 
Report  2:  Interactive  Graphics  Options 

Validation  Report:  Computer  Program  for  Design  and  Analysis  of 
Inverted-T  Retaining  Walls  and  Floodwalls  (TWDA) 

User’s  Guide:  Computer  Program  for  Design  and  Analysis  of 
Cast-in-Place  Tunnel  Linings  (NEWTUN) 

User’s  Guide:  Computer  Program  for  Optimum  Nonlinear  Dynamic 
Design  of  Reinforced  Concrete  Slabs  Under  Blast  Loading 
(CBARCS) 

User's  Guide:  Computer  Program  for  Design  or  Investigation  of 
Orthogonal  Culverts  (CORTCUL) 

User's  Guide:  Computer  Program  for  Three-Dimensional  Analysis 
of  Building  Systems  (CTABS80) 

Theoretical  Basis  for  CTABS80:  A  Computer  Program  for 
Three-Dimensional  Analysis  of  Building  Systems 

User's  Guide:  Computer  Program  for  Analysis  of  Beam-Column 
Structures  with  Nonlinear  Supports  (CBEAMC) 

User's  Guide:  Computer  Program  for  Bearing  Capacity  Analysis 
of  Shallow  Foundations  (CBEAR) 


Date 
Feb  1978 
Mar  1979 

Jan  1980 
Jan  1980 

Feb  1980 

Mar  1980 

Jun  1980 
Jun  1982 
Aug  1983 

Dec  1980 
Dec  1980 

Dec  1980 
Dec  1980 

Dec  1980 

Feb  1981 
Mar  1981 

Feb  1981 

Mar  1981 

Mar  1981 

Mar  1981 
Aug  1981 
Sep  1981 
Jun  1982 

Jun  1982 


(Continued) 
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(Concluded) 


Instruction  Report  K-83-1 
Instruction  Report  K-83-2 
Instruction  Report  K-83-5 

Technical  Report  K-83-1 
Technical  Report  K-83-3 

Technical  Report  K-83-4 
Instruction  Report  K-84-2 


Title 

User's  Guide  Computer  Program  With  Interactive  Graphics  for 
Analysis  of  Plane  Frame  Structures  (CFRAME) 

User's  Guide:  Computer  Program  for  Generation  of  Engineering 
Geometry  (SKETCH) 

User's  Guide:  Computer  Program  to  Calculate  Shear,  Moment, 
and  Thrust  (CSMT)  from  Stress  Results  of  a  Two-Dimensional 
Finite  Element  Analysis 

Basic  Pile  Group  Behavior 

Reference  Manual:  Computer  Graphics  Program  for  Generation  of 
Engineering  Geometry  (SKETCH) 

Case  Study  of  Six  Major  General-Purpose  Finite  Element  Programs 

User’s  Guide:  Computer  Program  for  Optimum  Dynamic  Design 
of  Nonlinear  Metal  Plates  Under  Bias*  Loading  (CSDOOR) 


Date 

Jan  1983 

Jun  1983 

Jul  1983 

Sep  1983 
Sep  1983 

Oct  1983 
Jan  1984 


